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We investigate the anisotropics in the cosmic microwave background in a class of models which 
possess a positive cosmic energy density but negative pressure, with a constant equation of state 
w = (p/p) < — 1. We calculate the temperature and polarization anisotropy spectra for both scalar 
and tensor perturbations by modifying the publicly available code CMBfast. For a constant initial 
curvature perturbation or tensor normalization, we have calculated the final anisotropy spectra as 
a function of the dark energy density and equation of state w and of the scalar and tensor spectral 
indices. This allows us to calculate the dependence of the tensor-to-scalar ratio on w in a model 
with phantom dark energy, which may be important for interpreting any future detection of long- 
wavelength gravitational waves. 



I. INTRODUCTION 

For many years jjj , arguments have been made in favor 
of the ACDM model because of its desirable property of 
having spatial flatness while accommodating a low mat- 
ter density; conditions which seem to be favored by cur- 
rent astrophysical observations. A few years ago, how- 
ever, an old alternative Q resurfaced and has been receiv- 
ing considerable attention ||. This alternative, which we 
will refer to as $CDM, borrows from the literature on in- 
flation and postulates a slowly rolling scalar field whose 
energy density has only recently become cosmologically 
relevant. This slowly rolling scalar field provides a con- 
stituent of the universe which is smooth below the hori- 
zon scaleQ and whose evolution can be made arbitrarily 
slow. This is usually parameterized by the equation of 
state, w = p/p, with the model becoming indistinguish- 
able from the cosmological constant in the limit w — > — 1. 
Gcncrically 1 > w > — 1 and w is typically a function of 
time. 

Of particular interest to us was a paper by Caldwell 
Pj which introduced a <E>CDM model with w < — 1. This 
form of dark energy, dubbed phantom energy also has the 
remarkable property of satisfying p^+p^ < 0, just as in 
the standard <I>CDM case. It was shown that this model 
is consistent with both recent observations and classical 
tests of cosmology, in some cases providing a better fit 
than the more familiar models with w > — 1. 

There are several known ways to achieve w < — 1 
e.g. Ref. ||, but Caldwell's model makes use of a non- 
canonical kinetic term. There is no obvious motivation 
from particle physics for considering such a radical ex- 



tension to the theory, and indeed the model has been 
criticized on these grounds ||. However we regard this 
as a useful toy model because it allows us to simply ex- 
tend the equation of state parameter space continuously 
below w = — 1 without the additional complications in- 
herent in modifications of the gravitational sector of the 
theory. (We shall discuss how our results depend on this 
detail later.) Observational constraints are often maxi- 
mally likely at the value w = —1 and for the purposes 
of maximum likelihood analysis it is better not to have 
the region of maximum likelihood bumping against the 
edge of the parameter space. Thus it is advantageous 
to extend the domain of the parameter w, even if one 
is primarily interested, on theoretical grounds, in values 
w > -1. 

One aspect which was not studied in |5|] was the effect 
of gravitational wave perturbations in the CMB, and in 
particular the dependence of the tensor-to-scalar ratio 
on the equation of state w. In this paper we study the 
tensor anisotropy spectrum in the phantom model. We 
show that the normalization of the tensor perturbations 
in the CMB is only weakly dependent on the equation 
of state. The behavior of the tensor-to-scalar ratio is 
therefore dominated by the w dependence of the scalar 
perturbations. 

We are interested in the tensor contribution to the 
anisotropy spectrum because it is, in principle, observ- 
able and directly connected to the energy scale of infla- 
tion H , a fundamental parameter of considerable impor- 
tance! Our calibration of the relation between the ampli- 
tude of the anisotropy in the CMB and this fundamental 
energy scale is clearly necessary to interpret any future 
observations of this signal. 



*The speed of sound for a scalar field is the speed of light. 
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II. THE MODEL 

Phcnomcnologically, the property of negative pressure 
and positive energy density can be achieved by consider- 
ing the (unorthodox) Lagrangian density 



(i) 



for the phantom component <fi. We are adopting the met- 
ric convention g^ v = diag (— a 2 , a 2 , a 2 , a 2 ) where a(t) is 
the cosmological scale factor. To lowest order the non- 
canonical negative kinetic term in the Lagrangian pro- 
duces the following expressions for the pressure and en- 
ergy density of this dark energy: 
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where over-dots represent derivatives with respect to con- 
formal time rj, and dt = adq. Negating the kinetic term 
also alters the equation of motion of the field from that 
of a field with a canonical Lagrangian by switching the 
sign of the dV/d<f> term. 
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Clearly the evolution of the scalar field is coupled into 
the evolution of the background Friedmann equations. 
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Here 7i = aH = a(a~ 1 da/dt), taking into account that 
our derivatives are with respect to conformal time 77. The 
effect of 4> on the evolution of the background cosmology 
can be implemented with only a slight modification to the 
CMBfast integration code M to include the new compo- 
nent's pressure and energy density. 

Once the expressions for p§ and are obtained, it is 
simple to invert them to obtain the potential, V in terms 
of the constant w and derivatives of the scalar field. In 
order to prevent tachyonic modes of the scalar field from 
developing, we address only models with constant w 0: 
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Following Q, by treating the dark energy as a perfect 
fluid, we can obtain an explicit expression for </>. 
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where p crlt = 3H 2 /8irGa 2 , and = p^/pcrit- Taking its 
derivatives we observe 
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Thus V and its derivatives will be functions of a and the 
parameters w, £1$, and po only. 

To calculate the anisotropy in the CMB we must in- 
clude the perturbations of <fi to first order. We work 
throughout in the synchronous gauge. Since we are solv- 
ing linearized equations it is useful to Fourier transform 
the perturbations so that we can solve the differential 
equations for each fc-mode independently. The equation 
of motion for a particular fc-mode is 
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where h is the trace of the metric perturbation in the 
synchronous gauge. The scalar field perturbation, being 
of spin-0, does not depend on the other metric perturba- 
tion directly. The 5(f>kS produce the fluctuations in the 
energy density and pressure of the (f> field: 
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These stress-energy perturbations must be included in 
the evolution of the metric perturbations for a self- 
consistent solution. The $CDM component contributes 
a source term to the right hand side of the Einstein 
equations. To make these calculations, we have modi- 
fied CMBfast by including the <&CDM contribution to 
the evolution of the background cosmology in the Fried- 
mann equations. We have also added the contribution of 
density and pressure perturbations of the <p field to the 
total curvature perturbation, using the results derived 
above. Even though the background pressure is nega- 
tive, we are certain that these perturbations are stable, 
because d 2 V/d<fr 2 is a negative definite quantity in this 
case. Examining Eq. (|l^) we see that the square of the 
effective mass will never become negative, and exponen- 
tially growing tachyonic modes will not occur. 

Note that 8<p cannot source gravitational wave modes 
because <j> is a scalar field, while a source of gravitational 
waves must have spin 2. Therefore, the only effect of 
(j) on the gravitational wave anisotropics is through the 
change in the evolution of the scale factor M . 



III. CMB ANISOTROPY SPECTRA 

Figs [l] and || display the scalar and tensor components 
of the temperature anisotropy for several values of the 
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equation of state parameter w. In each case we assumed 
that the cosmological parameters were h = .67, Vt m = 
0.3 = 1 — O^, fib = 0.04 and that the scalar spectral 
index n s = 1, and the tensor spectral index n t = 0. It 
is standard practice || to write the temperature field on 
the sky as a sum of spherical harmonics with coefficients 
ae m - The multipole moments, Ci, are then defined as 
Ce = \a,(, m \ 2 and the power spectrum, £(£ + l)Ce/(2n), is 
approximately the power per logarithmic interval in I. 

Since we are working to linear order in perturbation 
theory both the scalar and tensor spectra have one free 
overall normalization (ignoring for a moment any con- 
straints imposed by e.g. inflation). In our calculations 
we fixed the scalar normalization by requiring that the 
initial spatial curvature in the total matter gauge flf] ] 
£ = 1. For the tensor normalization^] we have chosen 
h = 1 initially. This choice makes it easy to interprate 
the dependence on w, and also to calculate the cosmol- 
ogy dependent fitting function presented in §IV. However 
with such a normalization the curves are very similar be- 
yond i ~ 100 for all values of w so for plotting purposes 
we have normalized the curves at £ — 10 in Figs [l] and 
||. In Fig. ||, the E-E and B-B (tensors only) polarization 
spectra are presented with Q = 1 or h = 1 normalization, 
which illustrates clearly that the height of the features is 
uniform with w. 



and the location of the peaks. The change in the height 
of the first peak traces to our decision to re-normalize 
the spectra to fit the COBE data. In fact the first two 
effects are manifestations of the same physical property: 
the suppression of the late Integrated Sachs- Wolfe (ISW) 
effect for lower w. In models such as $CDM, a period 
of accelerated expansion begins in the universe when the 
scalar field <f> becomes the dominant form of energy den- 
sity. For a fixed matter density Q m (measured relative to 
the critical density) the value of the parameter w deter- 
mines how early matter-(/> equality occurs, and therefore 
how long the universe has been experiencing the current 
period of accelerated expansion. When an accelerated 
expansion is occurring, potential wells which had formed 
earlier as a result of gravitational collapse begin to de- 
cay. If the well is decaying, photons which fall into the 
well gain more energy than they lose climbing out, which 
generates an anisotropy. This is the dominant effect at 
low I. For w < — 1 models, equality happens even later 
than in the case of the cosmological constant, while for 
w > —1, it happens earlier. For very negative values of 
w, there has been almost no acceleration in the recent 
history, so the potential wells have not decayed as dra- 
matically as they do in ACDM models. As a result, the 
late Integrated Sachs- Wolfe (ISW) effect has been almost 
entirely turned off B. 
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FIG. 1. The scalar CMB temperature anisotropy spectra, 
£(£+l)Ce/2(n) vs.£ for several values of w. The solid lines are 
w = —1.1, dashed are w = —3, dot-dash are w = —10, and 
dash-dash are w — —0.4. The spectra have been normalized 
at Cio to fit the COBE data. 

Consider first the anisotropy spectrum in Fig. [lj There 
are three features exhibiting w dependence: the height 
of the first peak, the height of the tail at low multipoles, 



^The tensor perturbation is gauge independent. 



FIG. 2. The tensor CMB temperature anisotropy spectra 
as a function of w. As in Fig. [l], the solid lines are w = —1.1, 
dashed are w — —3, dot-dash are w — —10, and dash-dash 
are w = —0.4. Notice that the curves are almost degenerate 
for the first 40 multipoles. There appears to be negligible 
dependence of the tensor normalization on the equation of 
state w. 

When normalized to £ = 1, the effect on the spectrum 
is to enhance the amplitude at low multipoles as w de- 
creases. Re-normalizing to fit the COBE data at £ = 10, 
the enhanced low-^ tail translates into a change in the 
peak height. For energy density perturbations, the de- 
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pendence of the spectrum normalization on the equation 
of state is principally due to the lack of late ISW effect 
as w decreases. 



subtends a smaller angle, shifting the anisotropy spec- 
tra to higher values of £. In the case of the polarization 
anisotropies, this is the only dependence on w (Fig. ||). 
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FIG. 3. The polarization anisotropy spectra. For the scalar 
modes, the E-E polarizations anisotropies are shown (top), 
and for the tensor modes the E-E and B-B polarization 
anisotropies are presented (middle and bottom). The scale 
is arbitrary, depending on our particular choice of normaliza- 
tion. 

The normalization of the tensor spectrum appears not 
to be affected by w in any significant way (Fig. ^|). This 
can be seen by noting that the peaks are all roughly at 
the same height, and the spectra between 2 < £ < 10 
agree closely (compare with Fig. [[]). The small w depen- 
dence that can be seen occurs for the lowest multipoles, 
corresponding to largest angular scales, and is again due 
to the ISW effect on the tensor spectrum at late times. 

For all the anisotropy spectra (both temperature and 
polarization anisotropies) there is a universal shift to the 
right for lower values of w. The reason for the shift is 
the same in all cases. The £ position of any feature in a 
spectrum depends upon the (comoving) angular diame- 
ter distance to the last scattering surface. This distance 
depends on w as follows: 



d\ ss = 



da 



where 



H 2 {a) = n r a 
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One notes that if w is more negative, this integral is 
larger, which implies that the last scattering surface is 
further away. A feature of a given physical size thus 



IV. DEPENDENCE OF T/S ON w AND n s 

For Qb — -04 and assuming a flat universe, we have 
explicitly calculated the evolution of the normalization 
of the spectrum from an initial curvature perturbation 
C = 1 and from h — 1. The normalization at the tenth 
multipole is equal to the value just after primordial infla- 
tion, multiplied by a cosmology dependent "transfer func- 
tion" , which we will call fs and fx- Using the method de- 
scribed in Jll| we have numerically determined this func- 
tion for the ranges —2.5 < w < —0.5, 0.6 < 0^ < 0.8, 
and —0.3 < nx < 0.3 and provide fitting functions be- 
low. The tensor spectral index nx, is related to the scalar 
spectral index nx — ns — 1 for power law models of pri- 
mordial inflation. We have deliberately left a w term 
in the evolution of the tensor normalization to illustrate 
how insignificant it it compared to the terms proportional 
to the spectral index and f^. A fit is presented below, 
where a = 0.6156,6 = 1.6, and c = —.2. The scalar trans- 
fer function has a mean accuracy of 2% over the range 
indicated above, and at worst is off by 10% near the edges 
of the parameter domains. This would correspond to a 
worst case 5% error in the temperature. Since the fit be- 
gins to go bad near the edges of the parameter space, we 
do not recommend extrapolating these functions further. 
The tensor function has a mean accuracy of .8%, with a 
worst case of 6% for a few points at the very edge of the 
fit. 

f s (w, nx,^) = -1.341 + 6mm a n b {~w) c + 1.936fi a 
+3.312n fc + 4.584(-w) c - 5.148fi a 7i b - 9.375Q a (~w) c 
-1.946n fc (-w) c + 5.172f7(-w) c - 2.692n(-w) c 

-.072n b w - .211(n- l) 2 (15) 



fx(w, n T ,^) = 0.358 + O.OOlw + 0.822n T - 0.046^ 
+0.337n^ + 0.029^! - 0A87n t Q 



(16) 



The tensor to scalar ratio is given by 
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(17) 



Fig. |i| shows the dependence of both the scalar and 
tensor normalizations at Cio on the equation of state w 
for an even broader domain than was used in the fitting 
formulae above. Observe that even over this huge range 
in the equation of state, the tensor normalization changes 
by less than 1%. 
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FIG. 4. The normalization of the tensor and scalar spectra 
at the 10th multipole for an initial perturbation £ = 1 and 
h = l. 



The longest wavelength modes contribute primarily to 
the very low-£ moments of the spectrum. A similar argu- 
ment can be made for the scalar modes, where the overall 
growth rate should be independent of the detailed model, 
but the longest wavelength modes will depend on any 
fluctuations in the "phantom" component through their 
impact on the evolution of the metric perturbations. Fi- 
nally, the angular shift with w of features in the spectrum 
is robust. 

We have quantified the dependence of the tensor to 
scalar ratio C[ /Cf on w, fi^, and tit- For the model 
of Ref . H , Eq. (O) allows us to relate the (potentially) 
observable C^/C^ to the energy scale of inflation, in the 
event that a gravitational wave signal is actually observed 
from the CMB. 
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V. DISCUSSION 

We have extended the work of Q on the $CDM 
model by calculating the tensor component of the CMB 
anisotropics. We have shown that the normalization of 
the tensor component is essentially independent of the 
equation of state w in this model. For a fixed value of 
the initial curvature perturbation after inflation £ = 1 
and a fixed initial gravity-wave amplitude, we have quan- 
tified the w, n<£ and tit dependence of both the tensor 
and scalar anisotropy normalizations at C\q for a broad 
domain of w. We displayed the anisotropy spectra for 
several values of w and discussed the features. 

Most of the results we present for w < — 1 are only 
weakly sensitive to the particular model chosen jl . While 
truly general statements are not possible, we can make 
several points about model dependence. The fact that 
our acceleration is driven by a field of spin zero implies 
there is no source of gravitational waves as <f> begins to 
dominate. This would remain true for other scalar (or 
vector) driven theories, while more complicated assump- 
tions (e.g. higher order gravity theories) may have ex- 
plicit source terms. However, we would expect only the 
longest wavelength gravitational waves to differ due to 
this in any significant way. This is because the accel- 
eration takes place very late in the history of the uni- 
verse for w < —1, so only the longest modes are affected. 
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